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SUMMARY 


Measurements of unsteady-flow turbulent boundary- layer characteris- 
tics have been made for artificially thickened boundary layers using a 
new optical technique which involves the bow-wave shape of bullets fired 
in a shock tube. Comparison of the unsteady-flow boundary -layer charac- 
teristics obtained with pyramid -shaped surface roughness is made with 
Schlichting's flat-plate steady-flow sand-roughness theory and with 
steady-flow measurements reported for sand roughness and for lateral 
square -bar roughness. 

Good agreement of the boundary-layer-thickness measurements with 
the steady -flow theory and with the steady-flow data was obtained. The 
measured velocity profiles were considerably less full than the steady- 
flow sand-roughness theory predicts but agreed well with steady-flow pro- 
files which were reported for the lateral bar roughness. It is shown 
that these differences in profiles do not lead to any large differences 
in the momentum thicknesses but do lead to much larger differences in the 
displacement thicknesses. 

Shock-wave-attenuation measurements were made in the shock tube for 
both smooth surfaces and surfaces of pyramid -shaped roughness of different 
heights and length. Comparison of the measurements with a simple atten- 
uation theory based on equivalent flat-plate steady-flow sand-roughness 
boundary- layer theory is made. 

The shock-wave-attenuation measurements agreed reasonably well with 
the theory for small values of the boundary-layer displacement thickness 
but the measurements were considerably lower than the theoretical values 
for larger values of the displacement thickness. 


INTRODUCTION 


In recent years there has been an increase in the use of the shock 
tube as a means for the study of fluid properties and fluid flows at the 
high temperatures for which this research tool is particularly adaptable. 
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This has naturally led to a corresponding increase in the study of the 
departures of the real fluid flow occurring in the shock tube from the 
flow indicated by any particular theory. These departures have, in 
general, been more marked at the higher temperatures and accompanying 
higher flow velocities and have been attributed to a number of causes 
including relaxation effects, changes of the fluid properties, and 
effects on the fluid-flow characteristics near the walls due to heat 
transfer and skin friction. 

The purpose of this investigation is to study the last effect, 
that of skin friction, by artificially increasing the magnitude of the 
friction through the use of surface roughness. In this manner, very 
thick boundary layers may be produced, and large shock-wave attenuations 
noted. Direct observations of the boundary -layer characteristics are 
made with a new optical technique which involves study of the bow-wave 
shape from a bullet fired in the shock tube. The shock-wave attenuations 
which are obtained with the use of surface roughness are of such magni- 
tude that the troublesome effects of imperfect diaphragm burst and 
instrumentation inaccuracies are reduced to second order. 

A review of the shock-tube literature shows that the aerodynamics 
of the real shock-tube flow in certain cases might be very different 
from that expected from the ideal theory. (See, for example, ref. 1 for 
treatment of the ideal shock- tube theory.) It was found in references 2 
to 4 that regions of flow in the shock tube which would otherwise have 
conditions constant in time exhibit continuous variation of pressure, 
density, and Mach number with time. Changes of the shock-wave strength 
and other flow parameters with distance traversed and with shock-tube 
pressure ratio, Reynolds number, and so forth are noted in refer- 
ences 1 , 5 , 6 , and 7 . Theoretical and experimental studies have been 
made to correlate experimentally obtained shock-tube flow measure- 
ments with viscous unsteady-flow theory. In these studies the wall 
boundary conditions are included by considering the flow as a quasi -one- 
dimensional flow, where the flow conditions are averaged across the 
channel at representative cross sections along the shock tube and the 
resulting values included in a one -dimensional theory. 

In references 1 and 5, one representative cross section was chosen, 
that being at the entropy discontinuity, and the mass-flow decrement at 
this cross section, as a result of the boundary-layer displacement thick- 
ness, was presumed to be equal to the mass-flow decrement at the shock 
wave at the same time. The displacement thickness was determined 
according to the Rayleigh laminar problem of the instantaneous accelera- 
tion to constant velocity of a flat plate in a fluid at rest. 

In an investigation (unpublished) at the Langley Aeronautical 
Laboratory in 1953 > the authors collaborated on a theory which, like 
references 1 and 5> used the mass -flow decrement due to the boundary-layer 
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displacement thickness at the discontinuity as the characteristic 
parameter. This mass-flow decrement, however, was assumed to equal the 
shock-wave mass -flow decrement which existed at a later time, this time 
being equal to the time required for a sonic wave to travel from the 
discontinuity to the shock. The displacement thickness was determined 
from flat-plate steady-flow turbulent empirical data at a flat-plate 
length corresponding to the flow distance from the diaphragm to the 
discontinuity. 

In reference 2, a theory was developed in which the entire hot-gas 
region between the entropy discontinuity and shock wave was represented 
as a quasi -steady flow by the choice of a coordinate system fixed to the 
shock. The Blasius laminar -boundary -layer solution was applied with the 
wall velocity nonzero, and the changes of momentum and energy due to 
skin friction and heat transfer were thus obtained. It was assumed that 
these effects resulted in the generation of waves which could overtake 
the shock and reduce its strength. 

An inclusive treatment of the shock-tube flow for the viscous, heat- 
conducting case has been done in reference 4 in which the complete flow 
in the shock tube is considered. Numerical evaluation of this theory is 
obtained by assuming that, at any point on the shock-tube wall, the drag 
and heat transfer are equal to the flat-plate steady-flow values at a 
distance from the leading edge of the plate equal to the shock-tube-f low 
distance at that point. The flow quantities in the boundary layer are 
averaged into the free -stream flow; by integration of the unsteady-flow- 
characteristics differential equations, including some linearizations, 
the flow parameters at any point or at any time in the shock- tube flow 
may be found. This method is such that any skin-friction curve - either 
laminar, turbulent, or combinations thereof - may be used. The first- 
order effects at the entropy discontinuity are also included. The 
treatment of the expansion fan timewise, however, was as a discontinuous 
front, although, for the pressure ratios covered by the experimental phase 
of the work, this treatment produced no large discrepancies (that is, 

1.94 to shock-wave pressure ratios). 

In the aforementioned theoretical treatments the correlation of 
experimental and theoretical results was, in some cases, not particularly 
good and, in the other cases, only a rather narrow range of shock 
strengths was covered experimentally. In all cases, the boundary -layer 
displacement thickness was assumed to be small in comparison with the 
shock-tube cross section so that effects of strearnwise pressure gradient 
on the boundary-layer development could be neglected along with any 
effects of area change on the shock-tube flow. 

In references 2 and 8, the correct treatment of the boundary layer 
for unsteady laminar boundary layers was used. It was shown that, for 
the laminar case, equivalent steady-flow flat -plate treatment would not 
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yield correct results nor would the Rayleigh treatment in the case of 
finite -strength shocks. There have been, however, no solutions for the 
unsteady turbulent-boundary-layer case; thus, the validity of an equiv- 
alent flat-plate representation is not known. 

It is felt that, by the present investigation, information may be 
obtained as to the validity of assumptions such as the equivalence of 
unsteady boundary layers to steady-flow flat-plate representations and 
the effect of thick boundary layers on the shock-tube flow. It is 
recognized, of course, that large surface roughness disallows the exist- 
ence of a laminar boundary layer and, in fact, the turbulent skin-friction 
dependency on Reynolds number is absent under most of the conditions 
encountered in the investigation. It is felt, however, that some funda- 
mental information as to unsteady-flow boundary -layer development may be 
obtained, despite this limitation, by comparing the resulting information 
with steady-flow flat-plate theory with these same limitations. 
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SYMBOLS 

velocity of sound 
cross-sectional area 
shock-tube width 
local skin-friction coefficient 
total skin-friction coefficient 
shock-tube perimeter 

bow -wave parameter, 

shock-tube height 
bow -wave incidence angle 
roughness size 

height above surface of a hypothetical wall which would give 
a shock-tube area reduction equal to average reduction due 
to roughness 

distance shock has traveled from leading edge of surface to 
midpoint of velocity -measurement section 
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fluid -flow Mach number 

reciprocal of velocity exponent in boundary layer. 



absolute pressure 

ux 

Reynolds number based on flow length, R = — 

■qJj 

Reynolds number based on sand roughness size, — - 

v 

time 

absolute temperature 
fluid-flow velocity 

shearing-stress velocity in boundary layer 

shock-wave velocity 

fluid-flow distance 

distance from wall normal to flow 


ratio of specific heats 
boundary-layer thickness 

incompressible boundary-layer displacement thickness (eq. (25)) 
incompressible boundary-layer momentum thickness (eq. (14)) 
roughness equivalence factor (eq. (l8)) 

Mach angle 
kinematic viscosity 
density 


wall -temperature coefficient. 
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Subscripts: 

0 high-pressure side of shock tube before burst 

1 low-pressure side of shock tube before burst 

2 region between shock-wave and entropy discontinuity 

00 free stream outside of boundary layer (as used herein is always 

that of region 2) 

a relative to fluid 

B bullet 

n normal to bow wave 

rw rough wall 

s sand roughness form 

saw saw profile roughness form 

sw smooth wall 

t total adiabatic 

w wall 

a immediately behind bow wave 


THEORY 

Bullet Theory 

The method usually adopted in steady flows for the determination 
of boundary -layer velocity profiles is the survey of this region by 
small pitot-static probes. This method is not adaptable, however, 
to unsteady flows since rapid-response gages of sufficiently small size 
are at present not obtainable. The use of an optical technique is 
desirable in general since the flow is not disturbed and photographs 
may be obtained at prescribed times. An interferometer could be used 
to determine density profiles and (for the case of constant static pres- 
sure) temperature profiles provided that refraction effects near the wall 
are sufficiently small to allow analysis of the records. In order to 
obtain velocity profiles, however, this method would depend directly on 
knowledge of the temperature -velocity relation in the boundary layer. 

The optical technique which is here set forth using bullets also requires 
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this knowledge of the temperature-velocity relation, but the results are 
not directly dependent upon the relation. The use of the bullet method 
also requires knowledge of the bullet bow-wave strength in the boundary 
layer, but it will be shown that lack of this knowledge may not be 
serious in many cases. 

If a bullet is moving supersonically with respect to a fluid, the 
disturbance, or shock wave, from the nose propagates out into the fluid 
and the path, or shape of this bow wave, is a history of the strength of 
the bow wave and of the fluid velocity and temperature (that is, fluid 
Mach number) through which it has propagated. If the wave should, at 
any point, encounter conditions tending to make the wave speed become 
subsonic with respect to the fluid, the shape is dependent also on time 
and eventually the wave disappears. In order, therefore, that interpre- 
tation of the bow-wave shape might be made in terms of velocity and 
temperature changes, the restriction of supersonic wave speed at all 
points is imposed and the change in wave strength must be known. 

Consider a bullet which is fired into a fluid in a direction opposite 
to the fluid flow (see following sketch) ; 
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The Mach number of the bullet with respect to the fluid is 



( 1 ) 


If streamwise variations in the fluid flow may be considered small in 
comparison with changes normal to the fluid flow, the Mach number normal 
to the wave is 


(it should be here noted that the wave incidence angle is a measure of 
the ratio of normal Mach number (wave strength) to the Mach number of 
the wave relative to the fluid only within certain limitations in unsteady 
flow, these limitations being that no gradients of u or a in the 
direction of flow be present within the region of interest. This region 
of interest (see preceding sketch) is that region bounded by the wave 
on the one side and by the propagation path of the lowest point of interest 
of the wave on the other side. If these conditions are respected, the 
incidence angle will then be a measure of the aforementioned ratio in the 
presence of gradients of u and a normal to the flow. ) 

If it can be assumed that the wave is very weak - that is, a Mach 
wave where = 1.0 - or if the wave strength is known, Mn is deter - 


If the wave strength is not known in the region of interest, some 
assumed variation of wave strength from a known point outside the region 
must be used. For example, the wave-strength variation in the free stream 
is known from determination of the wave angle, up, Uco, and a<», and 
this variation may be assumed to extrapolate into the region of interest. 

The bullet velocity ug is measured experimentally and the shock 
angle i is measured from schlieren photographs; thus, there are two 
remaining unknowns u and a and it is necessary to find an additional 
relationship. The relationship between the velocity and temperature in 
the boundary layer may be determined for particular cases. 


Mn = Ma sin i ^ 1.0 


( 2 ) 


mined since 



(3) 
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The simplest assumed relation between velocity and temperature in 
the boundary layer - for the case of no heat transfer at the wall and 
a Prandtl number of 1 - is used so that, according to the energy equa- 
tion, the relation between velocity and temperature in the hot-gas 
region may be written 





oo 


For the purpose of finding the influence of this assumption of no 
wall heat transfer on the results, the velocity-temperature relation for 
the turbulent case of a Prandtl number of 1, including heat transfer 
from the wall, will also be used for comparison. According to Crocco 
(ref. 9) > this relationship may be expressed as 

_T_ ^w _ T 2,t ~ T 2 /u\ 2 T 2,t ~ T w _u_ 

T2 T2 T2 \U2/ T2 U2 

which for the boundary condition = ^-w identical to equation ( 4 ) . 

For this unsteady case, however, the condition T w = is used. Equa- 
tion (5) is for the case of constant wall temperature and no streamwise 
pressure gradients. 


It is now possible to find the velocity at any point in the flow in 
terms of the bow -wave angle at that point if the wave strength is also 
known at that point. By combining equations (l) , (2), and ( 4 ), the 
solution for the velocity may be expressed nondimens ionally for the 
case in which there is no heat transfer from the wall as 


u 

^2 



2 (F + l) 


( 6 ) 


F = 


2 /sin i\^ 

T^\KT) 


where 
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By combining equations (l), (2), and ( 5 ), the solution for the 
velocity may similarly be expressed for the turbulent case of heat 
transfer from the wall only as 


u 

u 2 



2(F + 1) 


(7) 


where 



For the case in which there is no heat transfer from the wall, a value 
of -1 for 0 w makes equation (7) again identical to equation (6). 


Boundary-Layer -Growth Theory 

S mooth plate . - For the purpose of determining a theoretical boundary - 
ayer growth with which experimental measurements may be compared, as 
well as for later use in a simple attenuation theory, the equivalent 
steady-flow concept without heat transfer or pressure drop will be used. 
That is, flat-plate steady-flow turbulent-boundary-layer theory will be 
fl ° V length x = ut in th e shock tube is equivalent to 


If the initial assumption is made that the flow is turbulent from 
the leading edge of a plate and that the velocity distribution remains 
constant aiong the plate, or more specifically that 9/6 is constant, 
the boundary-layer growth along the plate may be readily determined in 
terms of the local skin-friction coefficient. (See, for example 
refs. 10 and 11.) ^ ’ 


d S(x) = 5 | c f dx 


( 8 ) 


If the velocity distribution in the boundary layer is of the power 

form 


u (y) 

Uco 


i/n 


y 

6 


( 9 ) 
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the local skin-friction coefficient for a smooth plate for Reynolds 
numbers in the range 5 X 10^ < R < 10 ? is found from reference 10 to be 


c f = 0.059R" 2/ ^ 5+N ) 


( 10 ) 


where N = 7- 

By substituting equation (10) into equation (8) and integrating 
from 0 to x, the boundary -layer growth along the smooth plate is 
found in terms of the total skin-friction coefficient: 


S(x) = I & xC f 
2 0 1 


where (see ref. 11) 


Cf = 


Cf 


3 + N 
1 + N 


( 11 ) 


( 12 ) 


A plot of the total skin-friction coefficient as a function of 
Reynolds number based on plate length according to equation (12) is 
given in figure 1. This curve is seen to approximate closely the smooth 
plate skin-friction value of reference 10 shown as the lower dashed line 

Equation (ll) may be expressed numerically for N = 7 in terms of 
the flow length x: 

5(x) = 0.38 xR -1 / 5 (13) 


since 


e n 

5 " (N + 1) (N + 2) 



( 14 ) 


with the assumption of low-velocity (incompressible) flow. 


Rough plate.- For fully developed roughness flow along a'flat plate 
that is, flow in which the roughness projects out of the laminar sub- 
layer, the steady-flow friction drag consists primarily of the drag of 
the individual roughness elements. In this case, the skin friction is 
independent of Reynolds number and for values of the characteristic 
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roughness number 



is given by 


C f = ^1.89 + 0 . 70 ^ logg ^ 


-2.5 


(15) 


(See ref. 10 for complete treatment of roughness flow over plates as 
transformed from Nikuradse's (ref. 12) measurements on pipes with sand 
roughness,) 

The velocity distribution for this case of fully developed roughness 
is assumed to be the same as that in Nikuradse's pipe and is given by 


£ " 2 -5 f s * 8.5 


(16) 


The above expression for total skin-friction coefficient (eq. (15)) 
is shown in figure 1 and the velocity distribution (eq. (l6)) is shown 
in figure 2. Note that the roughness skin friction in figure 1 is also 
shown as a function of flow Reynolds number for given values of the 
roughness Reynolds number Rks* (it is seen from fig. 2 that the 
Schlichting velocity distribution is very closely approximated over much 

of the range of -3- by an N = 5*5 power velocity profile except for 
K s 

the region close to the roughness.) 

Equation (l6) may also be expressed as 


u 

Uoo 




Uco 


2 ' 5 loE e | + 8 -5 TT 


h 

Uoo 


(17) 


Roughness other than sand may be used, however, with the same 
relations if an equivalent sand -roughness size is determined from 


k s = Ak 


(18) 


where k denotes the arbitrary roughness size and the factor A pre- 
scribes an equivalent sand size as used by Nikuradse to produce the same 
friction as that produced by k. (See ref. 13 for comparison of friction 
measurements of other roughness forms with sand. ) 
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The boundary-layer growth along the rough plate for fully developed 
sand roughness is found from substitution of equation (15) into equa- 
tion (ll) : 


5 


- 2.5 


( 19 ) 


Numerical evaluation of 


as a function of - — is now obtained 
k s 


8 u^ 

by first assuming a value of — in equation (17) and solving for — 

K s U-r 

U e 

at the edge of the boundary layer. If these values of — and — and 

Ur k s 

0 

equation (17) are used, — may be found from integration of equation (it) 
Equation (19) can then be plotted in the form of 


2 (— = ^L/i.89 + 0.704 loge " 2 ’ 5 


(k s 5 


0 


against t— and may be used to find and being known. This 

k s K s o K s 

function, which is the growth of the momentum thickness along the plate, 

is plotted in figure 3. 


It is recognized that equation (19) was obtained on the basis of 
a velocity distribution along the plate invariant with x and that 
this qualification does not hold for the case of roughness. (See 
eq. ( 17 )-) Since the change in velocity distribution along the plate 
has only a small influence in this case (see fig. 4 for influence of 
5 0 

on — according to the Schlichting theory) , equation (19) is justi- 
K s O 

fiable in view of some of the other assumptions which are made. 


Shock-Attenuation Theory 

A simple shock-attenuation theory will be developed for comparison 
with experimental results. It is similar to those of references 1 and 5 
in that the mass -flow decrement due to the boundary-layer displacement 
thickness at the entropy discontinuity will be presumed equal to the 
shock-mass-flow decrement at the same time. The theory is different, 
however, in that the boundary- layer displacement thickness will be 
based on the equivalent flat-plate steady-flow value for a plate of 
length equal to the shock- tube-flow length in the same manner as was 
done for the boundary -layer growth. 
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The reduction of shock-tube cross section at any location due to 
the boundary-layer displacement thickness is 


AA 5*(D - 45*) 
A A 


( 20 ) 


where D = 2b 4 - 2h is the shock- tube perimeter and A, the cross- 
sectional area. 

When 4b*«D, equation (20) reduces to 


AA( 5*«D) 

A A 


( 21 ) 


which is applicable when the displacement thickness is small in com- 
parison with the shock- tube cross-sectional dimensions. When this is 
not the case, the error in area reduction by using equation (21) would 
be 


AA( 8*«P) _ i = 1 

AA D . 

£&* “ 1 


( 22 ) 


If the free -stream conditions are assumed to be unaltered by the 
presence of the boundary layer, the fractional mass-flow reduction is 
then equal to the fractional area reduction: 


A(P 2 u 2 ) AA 

p 2 u 2 A 1 ' >) 

In order to convert the mass -flow reduction to shock-wave attenua- 
tion, it is assumed that the mass flow is reduced at the shock wave by 
the same fraction of the total as it is at the entropy discontinuity. 
(The wave mechanism by which this mass-flow reduction is accomplished 
is not considered, but the net effect would be that of expansion waves 
generated at the wall which would in turn overtake the shock wave.) 

The mass flow per unit area accelerated by the shock wave is nearly a 
linear function of the shock-pressure ratio (see fig. 5) j thus, the 
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shock-pressure attenuation is simply expressed for small values of 5* 

t>y 



Pi 


P2 

Pi 


1 



(24) 


For the smooth-wall incompressible turbulent boundary layer with a 
power-law velocity profile, the ratio of 6* to S is found from 


T = (u ”- u) d l = FTT 


(25) 


so that substitution of equations ( 25 ) and ( 13 ) into equation (24) for 
N = 7 yields 

. P2 


*2-1 

Pi 


0-58 D /_V| 1/5 
et - a 1 uox ) 


X 


( 26 ) 


and from the shock-tube x-t diagram (see fig. 6(a)) 


x , u 2 

x = T ! = vj 1 


(27) 


Equation (26) may now be rewritten as a shock-pressure attenuation 
coefficient : 


4 L 


A — 

Pi 


D P2 
Pi " 


Ur 


, 0.8 


k a l. 


" 0 T 2 oTB 
I ( Z sw) 


= 0.19 


R^0.2 /yqV 0.2 


( 28 ) 
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The right-hand side of equation ( 28 ) is a function only of shock- 
pressure ratio and initial temperature in the shock tube. 

For the rough -wall case the ratio of 5 * to 5 may be found by 
integration of equation ( 25 ) after substitution of equation ( 17 ) for the 

5 0 

assumed values of — , as was done in the case of — previously. 

k s o 

Substitution of equation (19) into equation (24) yields 


A p l = W 6D 
Po 2 6 9 A 

— - 1 
P 1 


^(1.89 + 0.704 loge tH 




(29) 


for the shock-pressure attenuation, and this may be expressed as a shock- 
pressure-attenuation coefficient by 


Pg 

Pl 


1 

^rw P2 

57 ' 


- 2 IT I TtI 1 - 8 ? + 0-T<* loge 57 -g- 


u 2 


-2.5 


l- rw\ 


8 9V 


(30) 


The right-hand side of equation (30) is a function only of the 


shock-pressure ratio and the ratio 


°rw 


Here again, the velocity 


distribution along the plate is not constant as the derivation requires, 

§* 0 

and the influence on — is more pronounced than it was on — (see 

o o 

fig. 4 ); thus, the error in attenuation will be more pronounced than was 

the error in boundary-layer thickness. 


APPARATUS AND PROCEDURE 


Studies of the boundary- layer growth and shock-wave attenuation were 

made in a stainless -steel shock tube 4 inches wide and 7 ^ inches high 

with a maximum length of 22 feet. A schematic arrangement of the 
experimental apparatus is shown in figure 7- The high-pressure side of 
the shock tube in this investigation was 4 feet long and pressures up 
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to 100 pounds per square inch were used. The low-pressure side con- 
sisted of interchangeable sections of 2-foot, 4-foot, and 8-foot lengths 
and of two sections each of 2-foot lengths, these latter sections con- 


taining 8- by l6 -inch plate -glass windows 1— inch thick. 


The sections 


were doweled for accurate alinement and connected with quick-acting clamps. 
The low-pressure side was at atmospheric pressure for shock pressure 
ratios up to about 2.5 and was lowered to about 0.1 atmosphere for the 
highest pressure ratios used. 


The smooth -wall attenuation study was made first after which the 
2-foot window sections were equipped with roughness on only the lower 
wall for the optical investigation using bullets. Finally, either the 
8-foot or 4-foot sections were alternately lined with roughness of 
three different grades on all four walls for the rough-wall attenuation 
studies. The roughness consisted of plastic sheets which were machined 
on the surface to form continuous rows of pyramid- shaped elements of 
square base and of height equal to one-half of the base dimension. The 
pyramids completely covered the surface and the rows were oriented at 
45° to the flow direction. (See fig. 7*) Three sizes of roughness were 
used but in each case the overall thickness of the plastic was such that 
the reduction of shock-tube cross-sectional area was the same for all 
three sizes. Smooth sheets of such thickness as to give this same area 
reduction were also used to line the 8-foot section in one series of 
tests to determine the effect of this area reduction on the shock-tube 
flow. All the sheets had 45° beveled side edges for fitting in the shock 
tube and were cemented tightly to the shock-tube walls. 

Shock velocities were measured by means of an electronic chronograph 
and two pickup stations were used to start and stop the chronograph. The 
pickups were located 18 inches apart in the 2-foot section and were 
commercial piezoelectric -type gages connected to the shock tube through 
slit-shaped openings on the walls of the shock tube. It was felt that 
this arrangement would minimize the voltage-rise time at shock passage 
since this rise time varies with the dimension of the opening in the flow 
direction. For the smooth-wall studies, the midway point of the gages 
was 15 feet from the diaphragm; for the rough -wall studies, the midway 
point was 1 foot downstream of the section containing the roughness. 

Records of pressure against time were obtained at the midpoint of the 
velocity gages in a few tests using capacitive -type flush-mounted pickups 
along with a cathode-ray oscilloscope and drum camera. 

For the bullet investigation of the boundary layer, 22-caliber "long 
rifle" and "Hornet" bullets were fired upstream (toward the high-pressure 
end) through a hole in the end plate of the low-pressure section. A 
paper diaphragm could be placed over the hole for evacuation purposes. 

The bullet speed could be approximately predetermined by the charge loading 
and, by means of pickups and delay circuits strategically placed, the 
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bullet speed could t>e electronically determined and synchronization of 
the optical light source with the bullet and shock-tube shock wave was 
possible. 

The optical components included a high-voltage spark-light source, 
thyratron controlled, and a two-mirror schlieren system using front - 
surfaced, l6-inch-diameter, parabolic mirrors. The spark duration was 
about 1 microsecond and there was sufficient light to expose a 5- by 
7-inch sheet of fast panchromatic film. 

The shock-tube diaphragm consisted of one or more sheets of exposed 
film so that the diaphragm would be stressed nearly to the bursting point 
and firing was accomplished with a solenoid-actuated trigger unit. 

The shock-tube temperature was determined prior to each run. 


RESULTS AND DISCUSSION 
Boundary- Layer Measurements 


Boundary -layer measurements were made with the bullet technique 
for two values of flow length over roughness and for two values of shock- 
tube pressure ratio (that is, flow velocity). The shock-tube configura- 
tions used for this investigation are shown in table I. It was necessary 
that proper synchronization of the shock tube, gun, and light source be 
accomplished so that the bullet would be within the "roughness -flow 
region" at the time of the photograph. That is, the photograph was taken 
at a time before any flow from the smooth region ahead of the roughness 
could reach the position under investigation (see fig. 6(b)) as well 
as before the shock reflected from the end plate could reach this posi- 
tion. The actual flow length appropriate to each photograph could then 
be determined from the x-t diagram using the reflected shock from 
the end plate as a reference point. The shock-tube flow was known from 
shock- velocity measurements. A photograph typical of those obtained 
using this method is shown in figure .8(a) . In order to determine the 
velocity profile from equation (6) or ( 7 ), it is necessary to know the 
Mach number normal to the bow wave or the bow-wave strength at each 
ordinate. The bow wave is nearly a Mach wave in the region of the 
boundary layer in these tests and its strength is closely determined 
outside the boundary layer from measurements of bow-wave incidence angle 
and the known bullet speed and free-stream conditions. (See eqs. (l), 
(2), and ( 3 ).) A plot of Mn as a function of distance from the bullet 
is then extrapolated into the region of the boundary layer. The coordi- 
nate system used in plotting the boundary -layer measurements is shown 
in figure 9 where the y = 0 point is taken at the ordinate corre- 
sponding to a smooth plane of equal shock- tube -area reduction. This 
system is also used in reference 13 . A plot of bow -wave incidence angle 
minus Mach angle near a rough wall as measured from a typical photograph 
(fig. 8(a)) is shown in figure 10. It can be seen from figure 10 that 
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the edge of the "boundary layer is readily determined in this case from 
the point at which the slope of the "bow wave changes more or less dis- 
continuously. It was found in many of the photographs, however, that 
the how wave was distorted as a result of interactions with one or more 
of the many weak reflections persisting in the field, these reflections 
having originated from the initial shock-wave movement over the roughness 
elements. (See fig. 8(h).) These distortions might appear anywhere 
along the "bullet how wave and, unfortunately, led to uncertainties in the 
analysis of the records. Selection of a more random type of roughness 
would prohahly have reduced this effect. 

In order that comparisons of the houndary-layer measurements might 
he made with the Schlichting theory, it is necessary to establish an 
equivalence of the roughness used to sand roughness. The theory outlined 
in Schlichting for a flat plate is based on measurements made by 
Nikuradse in pipes using sand roughness. For a roughness of different 
form and distribution from that of Nikuradse f s sand roughness, the 
equivalence of the two must he determined from friction measurements, 
as was done in reference 13 for many different types of roughness. Since 
the friction was not measured in this case, it was necessary to assume 
arbitrarily an equivalence criterion based on some physical concept even 
though there is no physical significance to the ratio A in equation (18). 
It was, therefore, assumed that surface roughnesses of arbitrary form and 
distribution would be equivalent on the basis of equal shock-tube-area 
reduction. From friction measurements given in reference 13, the equiva- 
lence of a saw profile roughness and Nikuradse ! s sand roughness was found 


The displacement thickness for the pyramid roughness used here is given 

by 


and, by using the assumed equivalence criterion, the pyramid equivalence 
is found by equating equations ( 31 ) and ( 32 ) to yield 



so that 



(31) 



(32) 


k - 8 k 
K s - Q K 


9 


(33) 
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With the equivalence factor A = the houndary-layer-thickness 

measurements obtained with the pyramid roughness are plotted in figure 11 
along with the Schlichting theoretical values from equation (19) • The 
agreement of the average of the experimental boundary-layer -thickness 
measurements with the theory is good although there is a noticeable 
scatter of the individual data points due to the aforementioned waves 
in the flow field. 


For comparison of these results with some steady-flow measurements 
obtained by using pitot-static survey probes, the results of Klebanoff 
and Diehl (ref. 14) for commercial sand roughness and of Moore (ref. 15 ) 
for lateral square bars are also plotted in figure 11. The equivalence 


of the Moore roughness form to Nikuradse sand roughness was taken as — 
^that is, A = and was determined by the same method as was used for 


the pyramids previously. The data of references 14 and 15 have been 
plotted by using the ordinate measured from y = 0 at the plane of 
equal average channel- area reduction. The agreement of the unsteady- 
flow boundary-layer-thickness measurements with those values of refer- 
ences 14 and 15 for steady flow is good but generally the points are 
lower than the theoretical values. It should be pointed out that a change 
in the roughness equivalence factor A would have only a relatively 
small effect on the comparison of the data with theory in figure 11. The 
values taken from reference 15 are for the roughness size closest to 
that used here. It was found in reference 15 that there was a slight 
dependence of the results on Reynolds number; this dependence, of course, 
is inconsistent with the sand-roughness theory. 


Velocity -Profile Measurements 


Boundary -layer velocity profiles over roughness were obtained from 
the bullet photographs and considerable scatter of the measurements was 
found as a result of the waves in the field, the influence being the 
most serious for the velocity profiles since few photographs were 
obtained without waves somewhere within the boundary layer. The average 
of the measurements taken at « 3 feet is plotted in figure 12 along 

with the theoretical distribution for this value of as determined 

k s 

from equation ( 17 ) . The experimental results were computed from equa- 
tion ( 6 ) . Also, for comparison, in figure 12 are plotted the steady- 
flow distributions found by Klebanoff and Diehl for commercial sand 
roughness and by Moore for lateral square bars at about the same values 

of — . The steady-flow data are plotted using the y = 0 point at the 
k s 

plane of equal channel-area reduction as was used for the reported data. 
It is seen in figure 12 that the boundary-layer velocity profile is 
somewhat less full for the pyramid and lateral-bar data (ref. 15 ) 
than for the sand-roughness case of reference 14, but that all 
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the experimental profiles are less full than the Schlichting theory 
predicts. Here again the influence of the roughness equivalence factor 
is very small. 

The most significant result appears to be that the non-sand roughness 
forms show a distinctly different profile from that of sand roughness 
although there is no apparent difference between the non-sand steady- 
flow and unsteady-flow profiles. In figure 13, the unsteady-flow data 
taken at x-^ w 1 foot are plotted in the same manner and compared with 
the Schlichting theory and the lateral -bar data of reference 15 for the 

same value of — . Here again, the two non-sand profiles are very 
s 

nearly the same although they are both less full than the Schlichting 
sand theory predicts. When the experimental unsteady-flow profiles were 
computed by using equation (7) , the results showed an even less full 
profile and gave a zero velocity at the peak of the roughness rather than 

at — k: therefore, these results were not used. 

3 

The profiles of figures 12 and 13 result in considerably larger 
experimental values of 5* than are predicted by the Schlichting theory 
as shown in figure 11. The experimental values shown result from inte- 
gration of the profiles in figures 12 and 13 according to equation ( 25 ) . 

5* 

The Schlichting theoretical curve for — is obtained by using figure 4. 

^ s 


Attenuation Measurements 

Shock-speed measurements including both smooth-wall and rough-wall 
studies were made by using the shock-tube configurations shown in table I. 
For each shock-speed measurement recorded, a shock-strength calculation 
was made from one -dimensional unsteady-flow shock theory, the measured 
room temperature being used and a value of 7 = 1.40 being assumed. For 
each shock-tube diaphragm burst the pressures on each side of the diaphragm 
were determined immediately prior to the rupture; from this shock-tube 

P 0 p p 

pressure ratio — , an initial ideal shock strength — was determined 

by using ideal shock-tube theory and 7 = 1.40. The difference between 
this initial shock strength and the shock strength determined from the 
velocity measurements is then considered as the shock loss due to the 
boundary layer only and is 
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and the shock-pressure attenuation is 




1 


For the purpose of determining a flow length x appropriate to each 
measurement for use in evaluation of the theory (see eqs. (26) and (29)), 
equation (27) is used where l is taken as the distance the shock has 
traveled from the leading edge of the surface under investigation to 
the midpoint of the velocity measurement section. (See table I.) The 

value of - is determined from the x-t diagram as the ratio of the 

distance a flow particle has traveled from the leading edge of the par- 
ticular surface to the distance that the shock has traveled from this 
leading edge during the same time. From ideal shock-tube theory this 
ratio is 


2 /P2 

x _ u 2 _ y - 1 \ p i 

l V]_ " 7 + 1 P2 

7 - 1 Px 

and is plotted in figure 14 for 7 = 1.4-0. 

Shock-pressure-attenuation coefficients were experimentally obtained 
for the last five configurations listed in table I from the left-hand 
side of equations (28) and (30) and the theoretical attenuation coeffi- 
cients were obtained from evaluation of the right-hand side of equa- 

tions ( 28 ) and (30) for the same — as before. The results for the 

k 

smooth-wall configuration are plotted in figure 15 and for the roughness 
configurations, in figures l 6 to 19 , all as a function of initial shock- 
wave strength. The solid arrows indicate the range of the root-mean- 
square deviation of the experimental data from the faired curve. The 
dashed arrows indicate the total range of deviation to be expected as a 
result of a 5 -microsecond error in time measurement in the velocity 
measurement section. The effect of poor diaphragm bursts on the data 
shown is not known but the influence of this factor on the smooth-wall 
data of figure 15 would be much more pronounced than on the rough-wall 
data, since the magnitude of the shock loss due to skin friction is of 
a different order. The same condition is seen to be true in the case 
of influence of errors in time measurement on the attenuation data. 
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The smooth-wall results shown in figure 15 are also representative 
of one series of tests not listed in table I in which smooth sheets 
8 feet long and having a thickness equal to the equivalent-area-reduction 
thickness of the roughness sheets were installed in the shock tube ahead 
of the measurement section. In view of this like result, it is felt that 
the 12 -percent area reduction due to the roughness does not appreciably 
affect the attenuation results; this result is probably due to the fact 
that the full shock-tube area is again attained at the measurement section. 

The number 4^ which appears in the shock-pressure-attenuation coefficient 

is the hydraulic diameter of the shock tube and the value appropriate to 
the results obtained when smooth or rough sheets were installed is, of 
course, smaller than the actual shock-tube hydraulic diameter. 

It can be seen in figure 15 that there is good agreement of the 
results with the simple theoretical concept used within the scatter of 
the data. However, the experimental scatter is such as to disallow a 
conclusive comparison. It can be seen in figures 1 6 to 19 that the 
agreement of the rough-wall results with the simple theory generally 
is good at the lower values of shock-pressure ratio but is increasingly 
poorer at the higher pressure ratios as the roughness height increases 
and as the roughness length increases. 


Pipe -Flow Effect 

Since the roughness covers all four walls of the shock tube and 
since the boundary layer does theoretically reach the center of the 
shock tube in a substantial part of the range of conditions covered, 
a pipe-flow effect is suggested; that is, as the ratio of 6 to the 

shock- tube half -width ^ approaches unity, the theoretical concepts 

which were used would be increasingly less applicable. The ratio -^7 

(or J increases with shock-pressure ratio, roughness size, and 

roughness length. In order to see whether any such pipe-flow effect is 
present, the data are shown in figure 20 as a ratio of experimental to 

2 6 * 

theoretical attenuation plotted as a function of — . The arrow 
indicates the abscissa at which 6 = ^b and to the right of which would 
be the pipe -flow regime. 


It can be seen from figure 20 that the data are generally less in 


agreement with the simple theory at the larger values of 
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The 


agreement of the data and theory for the high-pressure-ratio ends of the 
curves (right ends) appears to be very definitely a function of the 
displacement-thickness parameter. This statement might also be said of 
the low-pressure-ratio data (that to the left) except for that obtained 
with the largest roughness. It appears more likely, however, that there 


2b 
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is a pressure-ratio effect in addition to the effect of displacement 
thickness on the agreement of the results. There is a definite differ- 
ence in the nature of the pressure records for the high-pressure-ratio 
shocks and for the low values. In figure 21 is shown the pressure 
record obtained with configuration 7 for each of these cases. The 
pressure following the shock is seen to rise gradually with time in the 
case of the low-pressure ratio and to fall with time in the case of high 
pressure ratio. 


Qualification of Results 


In drawing any conclusions from the results which have been obtained 
from boundary -layer measurements and shock-attenuation measurements, 
certain points must be kept in mind. The boundary-layer thickness 6 
for the unsteady-flow pyramid roughness agreed well with both the steady- 
flow sand-roughness theory and the steady-flow measurements for sand 
roughness and lateral square-bar roughness. The velocity profiles for 
the unsteady-flow pyramid-roughness cases agreed well with the steady- 
flow lateral-bar cases when equation (6) was used but did not agree with 
the commercial sand-roughness case or the Schlichting sand theory. It 
cannot be concluded whether the agreement of the pyramid and lateral- 
bar data is due to the difference in roughness form being offset by the 
difference in the steady and unsteady flow or whether the two forms are 
equally different from sand and there is no difference due to the time 
nature of the flows. This question could best be settled, of course, by 
steady-flow boundary-layer measurements using pyramid roughness by methods 
similar to those used in reference 15 . It should be remembered when 
any comparison of the data shown is made with the theory that the 
Nikuradse profiles were obtained in pipes and represent streamwise 
pressure gradients whereas the other cases represent little or no pres- 
sure gradient. 


In any case, examination of equation ( 19 ) and figure 11 shows that 

q 

the effect of these differences in velocity profiles on the parameter — 

x 

and hence Cf (see eq. (ll)) is not large. 


The effect of velocity-profile differences on the boundary-layer 
displacement thickness (and hence on shock attenuation), however, is 
more pronounced, as can be seen from figure 11 and equation ( 29 ) . This 
effect is also seen in figure k where the values of the boundary-layer 


parameters 

o 


and — 
o 


that result from integration of the experi- 


mental profiles are shown. Therefore, the equivalence of an arbitrary 
roughness to sand rougjiness on the basis of equal friction is not suf- 
ficient for equivalence of the displacement thickness. It would be 
expected from this result that the attenuation theory which has been 
presented would predict values of attenuation which are too low. How- 
ever, quantitative refinement is not warranted in view of some of the 
other assumptions made in the simple theory. These assumptions include: 
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incompressible flow, constant velocity distribution along the surface, 
no wall heat transfer, a Prandtl number of 1, and use of boundary-layer 
characteristics at only the entropy discontinuity station to describe 
the attenuation. 


CONCLUDING REMARKS 


Measurements of unsteady-flow turbulent boundary -layer character- 
istics using the bullet technique and shock-wave attenuation in the shock 
tube with pyramid -shaped surface roughness and subsequent comparisons 
with steady-flow theory lead to the following conclusions: 

1. Good agreement of the boundary-layer-thickness measurements with 
flat-plate steady-flow sand-roughness theory and with steady-flow data 
using both sand roughness and lateral square bars was obtained. 

2. Measurements of the boundary -layer velocity profiles show a 
considerably less full profile than the steady-flow sand-roughness theo- 
retical profile but agreed well with steady-flow profiles for lateral 
square-bar roughness. It is shown that these differences in velocity 
profiles do not lead to large differences in momentum thickness but lead 
to much larger differences in the displacement thicknesses. 

3 . Shock-wave -attenuation measurements agreed reasonably well with 
a simple attenuation theory for small values of boundary-layer displace- 
ment thickness, but the measurements were considerably lower than theory 
at the higher values of displacement thickness. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 12, 1955- 
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TABLE I.- CONFIGURATION OF LOW-PRESSURE SIDE OF SHOCK TUBE (HIGH-PRESSURE 


SIDE - b FEET LONG IN ALL CASES) 


|~Air-to-air room temperaturej 







Distance 

from diaphragm, ft 



Number 

Configuration 

Roughness 

condition 

k, in. 

Smooth 

wall 

Rough 

wall 

Window 

section 

Shock 

velocity 

measure- 

ment 

section 

l 

1 

Bullet 

Roughness on 
bottom wall only 

0.25 

0 to 10 

10 to 16 

12 to 14 

8 to 10 

x rw ~ 

3 

2 

Bullet 

Roughness on 
bottom wall only 

.25 

0 to 14 

14 to 16 

14 to 16 

12 to 14 

x rw ** 

1 

3 

Smooth-wall 



0 to 14 



12 to 14 

13 


attenuation 






4 

Rough -wall 
attenuation 

Roughness on 
all four walls 

.0625 

0 to 8 
and 

8 to 12 


12 to 14 

5 






12 to 14 






5 

Rough-wall 

attenuation 

Roughness on 
all four walls 

.125 

0 to 12 
and 

12 to 16 


16 to 18 

5 






16 to 18 






6 

Rough -wall 
attenuation 

Roughness on 
all four walls 

.125 

0 to 4 
and 

4 to 12 


12 to 14 

9 






12 to 14 






7 

Rough-wall 

attenuation 

Roughness on 
all four walls 

.25 

0 to 4 
and 

4 to 12 


12 to 14 

9 






12 to 14 
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Figure 1.- Turbulent total-skin-friction coefficients for the smooth flat 
plate and the sand-roughened flat plate for incompressible flow according 
to Schlichting (ref. 10 ). Dashed region indicates the roughness flow 
which is not fully developed. 
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Figure 2. - Boundary-layer -velocity distribution for fully developed sand- 
roughness flow according to Schlichting (ref. 10 ). 
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Figure 3*- Theoretical growth of boundary -layer momentum thickness along 
a plate for fully developed sand-roughness flow. 
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Figure 4.- Comparison of theoretical and experimental dependency of the 
boundary -layer parameters on the boundary-layer thickness for fully 
developed roughness flow. 
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Figure 5« - Mass flow per unit area accelerated by a shock wave moving 
into a stationary fluid as a function of shock strength. 
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(a) Shock tube x - t diagram for low-pressure side showing the flow 
lengths appropriate to the smooth-wall and rough-wall flow. 



(b) Physical representation of boundary-layer growth along shock tube 

at time t sw . 


Figure 6.- Shock-tube flow representations. 
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Figure 7«- Schematic arrangement of shock tube and instrumentation 
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(a) Bullet bow -wave distortion at the rough wall. 



L-9l69ll 

(b) Reflections of the shock from the roughness elements. 

Figure 8. - Schlieren photographs typical of those obtained for the 

roughness conditions. 



Leading edge of roughness 


Figure 9 .- Application of the bullet technique to rough-wall boundary-layer 
measurements showing location of coordinates used. 
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Figure 10. - Measurements of the bullet bow -wave incidence angle minus 
Mach angle through the boundary layer on a rough wall. Ma = 1 *55« 
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Figure 11.- Comparison of experimental and theoretical boundary-layer 
thickness and displacement thickness along a rough plate for unsteady 
and steady flow. 
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Figure 12. - Comparison of experimental and theoretical boundary -layer 
velocity distribution for unsteady and steady roughness flow. 

8/k s ^ 10. 
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Figure 13.- Comparison of experimental and theoretical boundary-layer 
velocity distribution for unsteady and steady roughness flow. 

5/k s = 4. 
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Figure 14.- Theoretical ratio of fluid velocity to shock-wave velocity as 
a function of shock-wave strength for a shock wave moving into a sta- 
tionary fluid. 7 = 1.40. 
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Figure 15. - Comparison of theoretical and experimental shock-pressure 
attenuation coefficient as a function of shock strength for smooth- 
wall flow. Z sw = 13 feet. 
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Figure l6.- Comparison of theoretical and experimental shock-pressure 
attenuation coefficient as a function of shock strength for roughness 
flow# ^rwAs = 1>080; ^ 2 rw = 5 fset# 
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Figure 17*- Comparison of theoretical and experimental shock-pressure 
attenuation coefficient as a function of shock strength for roughness 
flow. Zrw/ks = 5^0? T-rw = 5 feet. 
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Figure 18. - Comparison of theoretical and experimental shock-pressure 
attenuation coefficient as a function of shock strength for roughness 
flow. Irv/kg = 972; l ^ ~ 9 feet. 
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Figure 19- - Comparison of theoretical and experimental shock-pressure 
attenuation coefficient as a function of shock strength for roughness 
flow# ^rw/^s = = 9 feet. 
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Figure 20.- Ratio of experimental to theoretical shock-wave attenuation 
as a function of the ratio of boundary- layer displacement thickness 
to shock-tube half -width. Arrow indicates abcissa at which boundary 
layer theoretically fills the shock tube. 
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(a) Weak shock. P 2 /P 1 = l-6l. 
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(b) Strong shock. P 2 /P 1 = 5*54. 

Figure 21.- Typical records of the variation of pressure with time of 
shocks after traveling over roughness. Configuration 7* 


NACA - Langley Field, Va. 












































